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COMPOSITE  FAILURES:  A  COMPARISON  OF  EXPERIMENTAL  TEST  RESULTS 
AND  COMPUTATIONAL  ANALYSIS  USING  XFEM 


INTRODUCTION 


Fiber-reinforced  polymers  (FRPs) — composites  consisting  of  a  thennoset  or 
thermoplastic  polymer  matrix  reinforced  by  carbon,  glass,  or  aramid  fibers — have  the  potential  to 
dramatically  influence  the  design  and  manufacturing  processes  of  innumerable  industries.  In 
applications  where  weight  is  a  premium  concern,  FRPs  present  unique  opportunities  to  improve 
existing  designs.  The  significance  of  FRP  composites  is  that  they  can  be  customized  in  the 
manufacturing  process  to  meet  the  needs  of  specific  applications.  Fiber  orientation  can  be 
tailored  to  provide  optimal  strength  and  toughness  for  certain  loading  conditions.  FRPs  are  of 
particular  interest  to  the  military  and  defense  industries  because  of  their  potential  use  in 
underwater  applications  where  corrosion-resistance,  weight,  and  implosion-resistance  are 
indispensable  design  factors. 

Because  of  the  widespread  availability  of  measurement  techniques,  experimental  testing 
of  composite  materials  has  largely  outpaced  the  computational  modeling  ability  of  such 
complicated  materials — thus,  relegating  the  design  of  composite  structures  to  a  build-test-build 
methodology.  Recent  advancements  in  computational  physics-based  modeling  (PBM),  however, 
are  changing  the  manufacture  of  FRPs:  comprehensive  modeling  techniques  can  significantly 
reduce  the  effort  required  to  build  and  test  future  composite  structures.  Previously  untested 
composite  arrangements  can  be  properly  analyzed  prior  to  any  physical  testing,  vastly  reducing 
cost  and  time  requirements.  Preliminary  computational  results  can  be  compared  to  experimental 
testing  results  to  validate  the  analytical  models. 

The  finite-element  method  (FEM)  has  been  a  popular  computational  analysis  technique 
for  modeling  mechanical  defonnations  and  failure  for  any  material,  including  FRPs.  Recent 
research  has  enabled  numerically  efficient  fracture  modeling  capabilities  using  the  extended 
finite-element  method  (XFEM).  Initiation  and  propagation  of  multiple  cracks  can  be  modeled 
simultaneously  without  mesh  refinement  or  successive  re-meshing  near  the  crack  tip.  Numerical 
computation  of  FRP  fracture  has  historically  consisted  of  imposing  failure  conditions  on  a 
homogenous  material  with  averaged  composite  properties  to  capture  one  of  the  common  failure 
methods,  such  as  delamination.  These  homogenized  materials,  however,  do  not  account  for  local 
effects  such  as  fiber/matrix  interfaces  or  fiber  orientation. 

This  report  documents  a  comparison  between  computational  analysis  and  experimental 
testing  of  unidirectional  (UD)  carbon  fiber/epoxy  composite  cylinders  under  lateral  compression 
between  two  flat  plates.  The  three-dimensional  (3-D)  computational  model  incorporates  fracture 
using  the  XFEM  capability  available  in  Abaqus  v6.13.  XFEM  is  used  to  model  crack  growth 
and  propagation  throughout  each  distinct  orientation-dependent  layer  of  the  composite,  with 
different  fracture  criteria  depending  on  fiber  orientation.  Experimental  and  computational  results 
for  crack  initiation  location  and  propagation  path  are  compared.  To  simulate  the  effect  of  surface 
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defects  in  experimental  specimens,  an  initial  crack  is  then  instituted  into  the  computational 
model.  To  further  demonstrate  the  influence  of  fracture  properties  on  analytical  results,  the 
fracture  parameter  inputs  of  the  computational  model  are  reduced  to  identify  the  effects  of  strain 
energy  release  rate  and  fracture  toughness  on  failure  response. 


FRACTURE  OF  COMPOSITES 


FRP  COMPOSITES 

FRP  composites  consist  of  a  thennoset  or  thennoplastic  polymer  matrix  reinforced  by 
carbon,  glass,  or  aramid  fibers.  Hundreds  of  thousands  of  individual  fibers,  whose  diameters  are 
5  to  10  pm,  are  bundled  together  to  fonn  a  tow.  Single-direction-fiber-tow  layers  can  be  (1) 
stacked  into  varying  orientation  plies  or  (2)  woven  together.  The  nature  of  manufacturing  FRPs 
provides  a  unique  ability  to  customize  the  FRP  for  specific  applications.  Fiber  and  matrix 
material  selection,  fiber/matrix  volume  fraction,  and  laminate  organization  are  just  a  few  of  the 
many  factors  associated  with  designing  and  manufacturing  a  desired  composite  material. 

Although  numerous  options  exist  when  the  use  and  design  of  FRP  composites  are  being 
contemplated,  the  exact  mechanical  nature  of  the  composite  has  not  been  fully  evaluated.  While 
FRPs  show  particularly  good  material  properties  in  the  fiber  axial  direction,  transverse  loads  on 
complex  structures  can  lead  to  failure  at  lower  stress  states.  These  issues  can  be  mitigated  by 
properly  designing  the  composite  layer  construction  to  provide  strength  in  strategic  component 
orientations. 

Typical  development  of  composite  prototypes  has  consisted  of  extensive  testing  of 
previously  built  structures.  For  example,  composite  cylinders  designed  for  underwater 
applications  are  first  subjected  to  rigorous  impact  and  implosion  testing  to  determine  their 
structural  efficacy.1,2  The  design  of  composite  structures  that  are  exposed  to  extreme  loading 
events  requires  that  attention  be  paid  to  material  failure  and  damage  absorption  characteristics.3 
Material  properties  can  be  derived  from  testing  of  circular  tubes,  and  the  analytical  solutions  can 
be  derived.4, 5  Cavallaro  researched  the  role  of  weave  styles  and  crimp  gradients  on  the  damage 
tolerance  of  woven  fabric  composites.6  The  majority  of  structural  research  as  it  pertains  to 
energy  absorption  capabilities  has  occurred  with  axial  compression  of  cylindrical  or  rectangular 
tubes,  where  local  wall  buckling  occurs.  Another  method  of  energy  absorption  is  through  lateral 
compression  of  thin-walled  cross-sectional  tubes.  Impact  energy  is  absorbed  as  strain  energy 
(elastic  and  plastic),  dissipative  energy,  and  fracture  energy.  To  properly  design  composites  to 
absorb  impact  energy  and  resist  damage  propagation,  fracture  mechanisms  in  composites  must 
be  well  understood. 
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FRACTURE  MECHANICS 


In  a  generic  material,  three  different  modes  of  fracture  can  occur  (see  figure  1).  A  review 
of  the  various  modes  of  fracture,  which  are  discussed  in  this  section,  is  provided  by  Prasad  et  al.7 

1 .  Mode  I  fracture  results  in  two  crack  faces  being  pulled  away  from  each  other. 

2.  Mode  II  fracture  occurs  because  of  in-plane  shear  loading,  where  the  two  fractured 
surfaces  slide  with  respect  to  each  other  along  the  crack  plane. 

3.  Mode  III  fracture,  out-of-plane  shear,  occurs  when  the  crack  faces  separate  out  of  the 
crack  plane. 


(a)  (b)  (c) 


Opening  In-Plane  Shear  Out-of-Plane  Shear 

Figure  1.  (a)  Mode  I,  (b)  Mode  II,  (c)  Mode  III  Fracture 


In  order  to  fully  understand  the  response  of  composite  materials  to  crack  initiation  and 
propagation,  fracture  mechanisms  must  be  researched  at  the  local  scale  through  controlled 
experiments.  Contrary  to  isotropic  materials  such  as  metals,  FRP  composites  possess  a 
discontinuous  fiber/matrix  interface,  where  complex  interactions  govern  material  failure.  As 
such,  FRP  composites  exhibit  localized  fracture  beyond  the  three  primary  fracture  modes  listed 
above.  Some  examples  of  composite  failure  are  (1)  fiber  fracture,  (2)  matrix  cracking,  (3) 
delamination,  or  (4)  fiber  pull-out.  The  most  common  methods  of  composite  failure  are  the 
failure  of  the  matrix  (matrix  cracking,  delamination)  or  the  failure  of  the  fiber/matrix  interface 
(fiber  pull-out).  Each  local  failure  phenomenon  relates  to  a  failure  in  one  of  the  three  primary 
modes  of  material  fracture. 

A  comprehensive  overview  of  the  mathematical  modeling  of  different  fracture 
phenomena  in  fiber  composites  is  provided  by  Mishnaevsky  and  Brondsted.8  In  addition,  Pinho 
et  al.9  developed  physically  based  criteria  for  calculating  multiple  modes  of  failure  in  FRPs. 
Depending  on  the  orientation  of  the  composite  layer  with  respect  to  the  advancing  crack  front, 
different  modes  of  fracture  occur  (see  figure  2).  For  cracks  that  are  parallel  to  the  fiber  direction, 
intralaminar  crack  growth  occurs  between  fibers  and  through  the  matrix.  If  the  crack  face  has 
formed  perpendicular  to  the  fiber  direction,  however,  crack  growth  requires  fiber  fracture. 
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Figure  2.  Crack  Growth  with  the  Crack  Face  Parallel  (Left)  and  Perpendicular  (Right) 
to  the  Fiber  Direction  (Perpendicular  cracks  require  fiber  fracture.) 


Fracture  toughness,  which  is  the  resistance  to  the  onset  and  growth  of  cracks,  is  a  key 
material  parameter  that  governs  the  ability  of  composites  to  absorb  damage.  This  toughness  can 
be  determined  by  comparing  the  total  strain  energy  release  rate  Gt  to  a  critical  material  fracture 
toughness  value  Gc ■  Gc  is  an  experimentally  determined  material  property  and  can  be  broken 
down  into  a  combination  of  the  individual  mode  fracture  toughness  values  Gic,  Gnc,  and  Gmc- 
Gt  is  a  summation  of  Gi,  Gu,  and  Gw,  which  correspond  to  energy  release  rates  due  to  mode  I,  II, 
and  III  fractures,  respectively.  When  GT  >  Gc,  the  crack  can  be  expected  to  propagate.  Using 
these  individual  strain  energy  release  rates,  an  overall  failure  criterion  can  be  established. 
Benzeggagh  and  Kenane  created  a  simple  mathematical  relationship  combining  the  various  mode 
strain  energy  release  rates  to  determine  Gc,  which  is  shown  in  equation  (l).10 

Gc  —  Gic  +  ( Gnc  +  G1UC  —  GIC )  (  " Gt 


Gt  -  G,  +  Gn  +  GIU.  (1) 

To  study  this  phenomenon,  Rebiere  and  Gamby  decomposed  the  strain  energy  release 
rates  for  different  modes  of  fracture  into  components  depending  on  the  type  of  crack  present; 
transverse  cracks,  longitudinal  cracks  or  delamination.11' 12  Pinho  et  al.  determined  the 
intralaminar  fracture  toughness  of  carbon/epoxy  composites  by  using  the  four-point  bend  and 
compact  tension  tests.13’ 14  Variations  in  mode  I  intralaminar  and  interlaminar  fracture  toughness 
have  been  observed.15, 16  Fiber  bridging  may  occur  in  interlaminar  failure  resulting  in  higher 
strain  energy  release  rates,  whereas  for  intralaminar  fracture  the  failure  path  navigates  between 
fibers. 


The  vast  majority  of  loading  conditions  on  composite  materials,  however,  involve 
mixed-mode  loading.  Even  when  global  single-mode  loading  is  applied  to  a  composite,  such  as 
double-cantilevered  beam  test,  localized  fiber/matrix  interactions  experience  mixed-mode 
loading.  Mixed-mode  bending  (MMB)  has  been  applied  globally  to  composites  to  ascertain  the 
mode  I  and  mode  II  delamination  strain  energy  release  rates.10, 17  The  present  authors  have 
investigated  localized  mixed-mode  fracture  for  single-  and  mixed-mode  loading  of  woven 
composites  as  well.18  For  curved  structures,  such  as  composite  tubes  under  lateral  compression, 
the  constituent  components  can  be  expected  to  experience  mixed-mode  loading. 
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EXPERIMENTAL  TEST 


TEST  SETUP 

For  composites,  fracture  mechanisms  can  have  significant  dependence  on  layer 
construction  and  orientation.  With  mixed-mode  loading  and  different  strain  energy  release  rates 
depending  on  fiber  orientation,  different  types  of  loading  can  produce  drastically  different  failure 
responses.  For  example,  composite  tubes  that  are  made  with  alternating  UD  composite  layers 
will  fail  differently  than  will  woven  composites  under  the  same  loading  condition.  Moreover, 
failure  modes  can  change  within  a  single-composite  construction,  depending  on  the  layer 
properties  impacted  by  the  crack  front.  To  demonstrate  the  different  failures,  lateral  compression 
was  performed  on  UD  carbon  fiber/epoxy  composite  tubes  composed  of  axial  0°  and 
circumferential  90°  direction  layers. 

A  cylinder  was  subjected  to  lateral  compression  between  two  flat  plates  until  primary 
failure  (see  figure  3).  The  test  was  performed  in  an  Instron  machine  with  a  1000-lb  load  cell  at  a 
rate  of  0.01  inch  per  minute  (0.25  mm/min).  The  loading  rate  was  selected  to  represent  a  static 
load  test  where  inertial  effects  can  be  neglected.  The  tested  cylinders,  produced  by  ACP 
Composites  (Livermore,  CA),  are  standard  modulus  wrapped  UD  carbon  fibers  in  an  epoxy 
matrix.  The  laminate  construction  consists  of  two  90°  hoop  layers  and  eight  0°  longitudinal 
layers. 


//////////////////// 


Figure  3.  Cylinder  Specimen  and  Experimental  Test  Setup  with 
Assumed  Failure  Location  at  the  3  o’clock  Position 
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Figure  4  illustrates  the  layup  construction  with  the  hoop-direction  layers  shaded  green. 
After  alternating  single-layer  fiber  orientations  at  the  inner  portion,  the  large,  outer  longitudinal 
layer  consists  of  seven  consecutive  layers  in  the  same  direction.  Four  specimens,  approximately 
3  inches  (76  mm)  in  length,  are  cut  from  the  same  cylinder.  The  inner  and  outer  diameters  are 
0.9  inch  (23  mm)  and  1  inch  (25  mm),  respectively.  For  a  cylinder  under  lateral  compression, 
maximum  principal  stress  concentrations  occur  at  the  vertical  inner  diameter  and  horizontal  outer 
diameter  locations.  With  longitudinal-oriented  fibers  on  the  outermost  layer  and  hoop-oriented 
fibers  on  the  innennost  layer,  fracture  is  expected  to  initiate  on  the  side  of  the  cylinders  at  the 
3  o’clock  or  9  o’clock  positions  (see  figure  3). 


Figure  4.  Partial  Cylinder  Specimen  Construction 
(Hoop  layers  are  shaded  green.) 
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TEST  RESULTS 


Failure  of  the  specimens  occurs  with  an  audible  crack  and  a  sudden  dropoff  of  force. 
Closer  inspection  reveals  that  the  failure  initiates  in  the  outer  layer  on  one  side  of  the  cylinder, 
matching  the  presumed  failure  location  described  in  figure  3.  After  failure  initiation,  the  crack 
propagates  radially  inward  through  the  matrix  until  it  interacts  with  the  reinforcing  fibers  of  the 
hoop  layer.  The  path  of  least  resistance  then  follows  along  the  interlaminar  interface  between  the 
hoop  and  longitudinal  layers  toward  the  3  o’clock  position.  The  failure  occurs  at  a  consistent 
position  near  the  3  o’clock  position  for  each  specimen,  where  the  maximum  tensile  stress  acts  on 
the  matrix  between  the  fibers.  Figure  5  shows  a  partial  view  of  a  failed  specimen  near  the  region 
of  interest.  A  series  of  enhanced  microscopic  images  of  the  crack  site  are  concatenated  in 
figure  6  to  highlight  the  crack  path.  Small  defects  and  delaminations  can  be  seen  in  the 
longitudinal  layers;  however,  it  is  unknown  if  those  defects  are  a  product  of  the  lateral 
compression  or  the  specimen  preparation  process. 


Figure  5.  Failure  of  Carbon  Fiber/Epoxy  Cylinder  Subjected  to  Lateral  Compression 


Crack  Initiation 

Intralaminar  Crack  Growth 


Delamination 


Figure  6.  Enhanced  View  of  Cylinder  Failure 
(After  mode  I  fracture  through  the  specimen,  the  crack  encounters 
a  reinforcing  hoop  layer  and  diverts  to  a  delamination.) 
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Figure  7  shows  the  load-versus-displacement  data  for  each  of  the  four  test  specimens. 
Each  specimen  dataset  is  adjusted  to  coincide  at  a  consistent  point  of  compression  where  a  linear 
load/displacement  curve  begins.  The  data  for  tubes  1,  3,  and  4  show  a  consistent  linear  slope; 
tube  2  experiences  gradual  cracks  before  failing  completely.  The  average  slope  of  force  versus 
displacement  for  tubes  1,  3,  and  4  is  5023  lb/in.;  whereas  the  slope  for  tube  2  is  4397  lb/in. 
Although  the  tests  were  stopped  at  this  point,  further  deformations  would  presumably  involve 
further  delamination  and  eventually  failure  of  the  hoop  fibers.  The  variation  in  failure  load  can 
be  attributed  to  the  presence  of  surface  defects  that  act  as  crack  initiation  sites. 


Displacement  (in) 


Figure  7.  Force  Versus  Displacement  for  Experimental  Results 


For  these  composite  tubes,  the  layer  construction  played  a  large  role  in  the  manner  of 
failure.  Longitudinal  layers  proved  to  be  weak  in  lateral  compression;  whereas  the  hoop  layers 
resisted  a  propagating  crack.  The  higher  failure  strength  and  fracture  toughness  of  the  fibers 
prevented  a  continued  radial  crack  path;  instead,  the  intralaminar  crack  was  diverted  to  a 
delamination  between  the  hoop  and  longitudinal  layers.  The  longitudinal  layers,  though  weak  in 
lateral  compression,  are  significantly  stronger  in  axial  bending.  This  bending  load  condition 
places  the  fibers  in  tension  and  compression,  giving  support  to  the  tube  structure.  Just  as  the 
longitudinal  fibers  are  weak  in  lateral  compression,  the  hoop  fibers  would  be  weak  in  axial 
bending.  One  method  of  improving  the  strength  and  damage  tolerance  of  these  tubes  under 
lateral  compression  is  to  add  a  hoop  layer  to  the  outer  skin  of  the  tube,  providing  circumferential 
strength  at  the  crack  initiation  point. 
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PHYSICS-BASED  MODELING 


FEM 


Because  of  the  widespread  availability  of  measurement  techniques,  experimental  testing 
of  composite  materials  has  outpaced  the  computational  modeling  ability  of  such  complicated 
materials.  Recent  advancements  in  computational  PBM,  such  as  FEM,  are  changing  the 
manufacture  of  FRPs.  Comprehensive  modeling  techniques  can  significantly  reduce  the  effort 
required  to  build  and  test  future  composite  structures,  reducing  cost  and  time  requirements. 

These  analytical  results  can  be  compared  to  experimental  testing  to  further  refine  existing 
designs  or  develop  novel  systems. 

The  majority  of  computational  modeling  has  focused  on  investigating  composites  in  a 
large  scale,  where  the  composite  is  assumed  to  be  a  homogenous  material  with  averaged 
mechanical  properties.  In  ordinary  static  analyses,  for  which  fracture  is  not  a  concern,  this 
assumption  may  be  sufficient.  For  dynamic  loading  and  other  analyses  that  focus  on  material 
failure,  localized  fracture  mechanisms  must  be  addressed.  Many  techniques  have  been  studied  in 
attempting  to  model  damage  in  composites.  Customized  damage  parameters  can  be  created  to 
initiate  failure  at  certain  thresholds.19, 20  Upon  reaching  the  threshold,  material  properties  are 
changed  to  reflect  the  loss  of  load-carrying  capacity  of  failed  elements.  Other  techniques 
incorporate  fracture  into  the  analyses.  Typically,  fracture  analyses  are  assessed  by  implementing 
an  initial  crack  into  the  component.'  Crack  growth  forms  along  a  predetermined  path. 

Numerous  studies  of  quasi-static  and  dynamic  impact  of  composite  tubes  have  been  investigated 
using  various  techniques.22’26 

Composite  materials,  however,  consist  of  multiple  distinct  materials  with  varying  failure 
mechanisms.  Drastically  different  material  responses  to  impact  occur  when  interactions  between 
the  fiber  and  matrix  are  accounted  for.  Segala  and  Cavallaro  have  investigated  the  failure  modes 
of  distinct  UD  fiber/matrix  composite  systems  subjected  to  ballistic  impacts,  blast  pressures,  and 
mechanical  shock  with  these  interactions  in  mind."  Even  still,  most  analyses  focus  on  consistent 
layer  properties  throughout  the  component  thickness,  which  neglect  the  orientation  dependence 
of  woven  and  UD  composites. 


XFEM 


While  most  of  the  research  on  the  fracture  of  composite  materials  has  used  conventional 
FEM  techniques,  the  XFEM  offers  unique  capabilities  in  modeling  fracture.  Though  standard 
finite-element  analyses  require  a  highly  refined  mesh  near  the  assumed  crack  location,  XFEM 
allows  cracks  to  fonn  and  propagate  along  a  natural  crack  path  without  prior  knowledge  of  its 
location;  moreover,  continual  re-meshing  near  the  crack  front  is  not  required.  Eliminating  these 
computationally  expensive  methods  makes  modeling  fracture  using  XFEM  markedly  easier. 
Widespread  research  has  been  conducted  implementing  XFEM  into  established  finite-element 
codes.'  Recently,  commercial  software  packages  such  as  Abaqus  have  included  an  XFEM 
capability.  With  these  new  abilities,  previously  impractical  simulations  involving  composite 
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structures  can  be  analyzed  quickly  and  efficiently.  In  turn,  XFEM  can  be  used  to  validate 
experimental  procedures  and  predict  future  test  results. 

This  analysis  utilizes  the  XFEM  available  in  Abaqus.  XFEM  has  been  widely  used  to 
model  fracture  in  composites  under  various  loading  conditions.  ’ "  '  An  initial  investigation  of 
the  efficacy  of  XFEM  was  previously  reported.  32  In  that  report,  XFEM  and  cohesive  zone 
modeling  (CZM)  were  combined  to  simulate  the  failure  mechanisms  of  distinct  fiber  and  matrix 
components  in  carbon  fiber/epoxy  tubes  under  hydrostatic  depth  pressure.  XFEM  was  shown  to 
adequately  identify  crack  initiation  locations  and  propagate  crack  growth  under  further  loading. 

Within  XFEM,  fracture  of  the  cylinder  is  governed  by  the  virtual  crack  closure  technique 
(VCCT)  under  the  assumption  of  linear  elastic  fracture  mechanics  (LEFM).  In  VCCT,  the 
energy  absorbed  through  material  fracture  is  assumed  to  be  equal  to  the  work  required  to  close 
the  crack  faces.  Krueger  documented  a  large  review  of  fracture  toughness  characterization  for 

T -5 

multiple  modes  of  fracture  using  the  VCCT.  XFEM  was  originally  developed  by  Belytschko 
and  Black34  and  is  based  on  the  partition-of-unity  method  (PUM)  developed  by  Melenk  and 
Babuska.  In  PUM,  additional  enrichment  functions  are  inserted  into  the  conventional  shape 
functions  to  account  for  discontinuities  across  the  crack  front.  This  continuity  allows  cracks  to 
form  and  propagate  without  prior  knowledge  of  the  crack  location  and  does  not  require  re¬ 
meshing  throughout  the  simulation  to  maintain  a  fine  mesh  near  the  crack  tip.  The  partition  of 
unity  is  further  described  in  detail  by  Belytschko  et  al.  The  applied  enrichment  functions 
introduce  additional  degrees  of  freedom  to  affected  elements,  allowing  the  element  to  undergo 
transformation  while  maintaining  continuity.  In  fracture  mechanics,  the  enrichment  functions 
consist  of  a  near-tip  asymptotic  function,  which  captures  the  stress  singularity  at  the  crack  tip, 
and  a  discontinuous  function  to  map  the  displacement  of  the  crack  faces  within  an  element.  The 
enriched  shape  function  is  given  in  equation  (2):36 

Uh(x)  =  £v,  N,  (x)U,  +  Y,rx;„  (x)  [«(/(*))  -  H  (/(*;))]  q° 

+  Eke sc  Nk O)  [¥«>(*)  -  V«>(xOk“’,U>.  (2) 

The  first  summation  tenn  includes  the  typical  finite-element  shape  function,  where  A7  are 
the  standard  shape  functions  and  uj  are  nodal  degrees  of  freedom.  The  second  summation  tenn 
adds  additional  enrichments  to  the  elements  in  the  set  Sh  encompassing  the  crack  tip.  H  is  the 
Heaviside  function  and  q{j  are  enrichment  coefficients.  The  third  summation  term  adds 
enrichments  to  the  elements  cut  by  the  crack,  excluding  the  elements  inherent  in  the  second  tenn 
and  are  included  in  the  set  Sc-  The  tenn  is  a  set  of  enrichment  functions  near  the  crack  tip. 
H  is  defined  in  equation  (3)  as: 

H{r>  =  ft  Otherwise'  <3) 
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Because  the  crack  only  affects  elements  through  which  it  passes,  only  the  connecting 
nodes  gain  the  additional  degrees  of  freedom.  Limiting  the  enrichment  function  to  nodes  of 
interest  maximizes  the  efficiency  of  instituting  the  additional  degrees  of  freedom.  Figure  8 
shows  a  schematic  of  how  the  nodes  are  assigned  around  a  given  crack. 


Figure  8.  Crack  Running  Through  a  Set  of  Elements  with  Enriched  Functions 


3-D  FINITE  ELEMENT  MODEL 


OVERVIEW 

In  order  to  properly  compare  failure  mechanisms  in  the  analysis  with  those  of  the 
experimental  specimens,  a  3-D  finite-element  model  was  created  in  Abaqus  to  represent  the  test 
cylinders.  Eight  node  brick  elements  (C3D8R)  were  utilized  to  construct  the  solid  model.  The 
cylinder  has  a  0.9-inch  (23  mm)  inner  diameter,  a  1. 0-inch  (25  mm)  outer  diameter,  and  is 
0.02  inch  (0.5  mm)  in  depth.  Geometrical  nonlinearities  are  accounted  for  in  the  static  analysis 
as  well.  To  reduce  the  number  of  elements  in  the  model,  a  half-symmetry  boundary  condition  is 
implemented.  The  half-symmetry  boundary  condition  is  imposed  on  the  vertical  midplane  due  to 
the  crack  location  at  the  3  o’clock  position.  A  long  cylinder  is  simulated  by  applying  symmetry 
boundary  conditions  on  each  z-direction  face.  Plane  strain  is  not  applicable  to  this  analysis 
because  of  the  different  modulus  of  elasticity  in  the  z-direction  for  each  layer.  To  match  the  test 
specimen  construction,  eight  layers  are  oriented  in  the  longitudinal  0°  direction  and  the 
remaining  two  layers  are  oriented  in  the  90°  hoop  direction.  The  cylinder  wall  thickness  is  0. 1 
inch  (2.5  mm),  with  each  layer  representing  a  0.01 -inch  (0.25  mm)  thickness.  Figure  9  shows  a 
side  view  of  the  model  and  load  condition.  The  hoop  layers  are  depicted  in  orange,  and  the 
longitudinal  layers  are  green. 


11 


Two  rigid  surfaces  represent  the  flat  plates  that  compress  the  test  specimen.  As  in  the 
experimental  procedure,  the  bottom  plate  is  held  fixed  while  the  top  plate  moves  at  a  prescribed 
vertical  displacement.  With  an  assumed  crack  location  that  is  based  on  the  experimental  results, 
only  a  subset  of  elements  in  the  model  must  include  the  XFEM  enrichment  functions. 


Figure  9.  Cylinder  Model  Construction  with  Hoop  (Orange)  and  Longitudinal  (Green) 
Oriented  Fibers;  XFEM-Enriched  Elements  Highlighted  in  Red 


MATERIAL  FRACTURE 

Each  layer  is  modeled  as  a  homogenized  material  using  a  linear  elastic  orthotropic 
constitutive  behavior.  The  fibers  are  oriented  in  the  local  1 -direction;  whereas  the  fibers  vary 
with  respect  to  the  global  1 -direction  depending  on  the  layer  orientation.  The  hoop  layers  are 
oriented  so  that  the  fibers  provide  support  to  the  tube  circumferentially.  This  direction 
corresponds  to  the  global  ^-direction  of  the  cylindrical  coordinate  system.  The  longitudinal 
fibers  provide  strength  in  the  global  z-direction,  which  is  not  subjected  to  significant  deformation 
in  this  analysis.  Table  1  shows  the  local  material  properties  used  for  both  the  hoop  and 
longitudinal  layers.  Although  the  detailed  material  properties  are  not  available  for  the  tested 
specimens  from  ACP  Composites,  the  applied  values  correlate  to  typical  carbon/epoxy  systems. 
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Table  1.  Constant  Material  Properties  for  the  Composite  Layers 


En 

E22 

E33 

G12 

G13 

G23 

ksi 

ksi 

ksi 

ksi 

ksi 

ksi  V12  V13  V23 

(GPa) 

(GPa) 

(GPa) 

(GPa) 

(GPa) 

(GPa) 

19000 

1300 

1300 

750 

750 

750 

a:  0.33  0.33  0.25 

(p./j 

(131) 

(9) 

(9) 

(5.2) 

(5.2) 

Damage  of  the  layers  is  controlled  by  a  maximum  principal  stress  condition.  The 
material  undergoes  fracture  above  the  prescribed  maximum  principal  stress  ops.  For  the 
hoop-direction  layers,  it  is  assumed  that  mode  I  crack  growth  would  result  in  fiber  fracture.  For 
the  longitudinal  layers,  the  crack  would  propagate  through  the  matrix  in  between  the  fibers.  As 
such,  the  hoop  layers  are  significantly  more  resistant  to  fracture  than  are  the  longitudinal  layers. 

The  Benzeggagh-Kenane  criterion  is  used  to  calculate  Gc  for  each  layer.  Once  the  total 
strain  energy  rate  Gt  exceeds  Gc,  fracture  is  assumed  to  propagate.  Table  2  details  the  strain 
energy  release  rates  and  principal  failure  stresses  for  the  hoop  and  longitudinal  layers,  which  are 
taken  from  literature.4' 5 


Table  2.  Fracture  Toughness  Values  and  Failure  Stresses 
for  the  Hoop  and  Longitudinal  Layers 


Property 

Gic  Guc  Gmc  GPs 

lb/in.  lb/in.  lb/in.  ksi 

(kN/m)  (kN/m)  (kN/m)  (MPa) 

Hoop 

Longitudinal 

140  140  1  15 

(24.5)  (24.5)  (0.175)  (103) 

1  1  1  350 

(0.175)  (0.175)  (0.175)  (2413) 

Using  XFEM  enables  the  modeling  of  fracture  for  both  the  hoop  and  longitudinal  layers 
without  modifying  the  parameters  of  the  initial  model.  Because  the  fibers  are  significantly 
stronger  than  the  matrix  and  have  a  higher  fracture  toughness,  the  longitudinal  layers  are 
expected  to  fail  before  the  hoop  layers. 


RESULTS 

Before  fracturing,  the  cylinder  experiences  the  maximum  principal  stress  in  the  innermost 
hoop  layer  near  the  plates  (see  figure  10).  This  location  is  a  combination  of  tensile  strain  as  the 
cylinder  bends  to  a  flatter  profile  and  a  high  modulus  of  elasticity  for  the  hoop-oriented  fibers. 
The  tensile  state  of  stress  is  similar  to  that  for  a  curved  beam  in  bending.  The  stress  at  this 
location,  however,  does  not  reach  the  considerably  higher  principal  failure  stress  of  the  fibers, 
and  fracture  is  not  achieved. 
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S,  Max.  Principal 
(Avg:  75%) 


+  1.228e+02 
+  1.124e+02 
+  1.019e+02 
+  9.149e+01 
+  8.105e+01 
+  7.061e+01 
+  6.017e+01 
+4.973e+01 
+  3.929e+01 
+  2.8B6e+01 
+  1.842e+01 
+  7.978e+00 
-2.460e+00 


Figure  10.  Maximum  Principal  Stress  Immediately  Before  Fracture  Initiation  (ksi) 


The  other  location  of  high  stress  is  at  the  outermost  layer  of  the  cylinder  at  the  3  o’clock 
position.  At  the  location  of  interest,  the  highest  principal  stress  occurs  in  the  outer  (longitudinal) 
layer.  Because  of  the  lower  modulus  of  elasticity  in  the  circumferential  direction,  the 
longitudinal  layers  experience  lower  stress  than  do  the  reinforcing  hoop  layers  for  the  resultant 
bending  strain.  At  the  region  of  interest,  however,  the  inner  hoop  layers  are  in  compression  and 
the  outer  layers  are  in  tension.  Figure  1 1  shows  the  maximum  principal  stress  contour  for  the 
cylinder  immediately  prior  to  fracture  initiation  near  the  3  o’clock  position. 


S,  Max.  Principal 
(Avg:  75%) 


H 


E 


+  1.500e+01 
+  1.354e+01 
+  1.208e+01 
+  1.063e+01 
+  9.167e+00 
+  7.708e+00 
+  6. 2506+00 
+  4.792e+00 
+  3.333e+00 
+  1.875e+00 
+  4.167e-01 
-1.042e+00 
-2.500e+00 


Figure  11.  Maximum  Principal  Stress  Immediately  Before  Fracture 
Initiation  at  Region  of  Interest  (ksi) 
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The  outermost  layer  of  the  matrix  reaches  its  failure  threshold  ops  before  the  vertical  inner 
hoop  layer  reaches  its  higher  failure  threshold,  inducing  crack  initiation  on  the  outer  edge.  The 
crack  is  first  directed  inward  and,  with  further  loading,  propagates  radially  inward  as  shown  in 
figure  12.  No  additional  cracks  are  formed  in  the  longitudinal  layers,  and  the  maximum 
principal  stress  in  the  hoop  layers  remains  below  ops.  At  a  displacement  of  0.125  inch 
(3.18  mm),  matching  the  failure  displacement  in  the  experimental  test,  the  crack  has  propagated 
only  partially  through  the  cylinder. 


Figure  12.  Partially  Propagated  Crack 


EFFECTS  OF  FRACTURE  PARAMETERS 

For  a  plate  displacement  that  causes  a  large  crack  and  delamination  to  form  in  the  test 
specimens,  the  computational  analysis  exhibits  only  partial  crack  propagation.  One  possible 
explanation  for  this  discrepancy  is  an  overestimation  of  the  strain  energy  release  rate  values.  A 
drawback  of  implementing  a  half  symmetry  boundary  condition  is  that  it  assumes  that  a  crack 
forms  on  both  sides  of  the  cylinder — which  contradicts  the  experimental  results.  The  additional 
strain  energy  introduced  by  modeling  a  full  cylinder  under  lateral  compression  could  lead  to 
further  strain  energy  release  via  crack  propagation.  Exploring  this  possible  energy  deficiency 
requires  reducing  the  strain  energy  release  rates  proportionally  in  both  the  hoop  and  longitudinal 
materials.  Varying  the  critical  strain  energy  release  rate  values  provides  a  method  for 
investigating  crack  growth  beyond  what  the  initial  model  parameters  can  exhibit.  At  a  lower  Gc, 
the  crack  requires  less  energy  to  grow.  Another  theory  for  the  decreased  crack  growth  is  the  lack 
of  an  initial  crack  site.  Realistically,  composite  materials  exhibit  numerous  defects  and  surface 
deformations  that  are  eligible  sites  for  crack  initiation.  In  this  finite-element  model,  the  energy 
that  is  needed  to  fonn  a  crack  must  first  be  achieved  before  crack  propagation  can  commence. 

To  facilitate  the  crack  growth,  an  initial  crack  is  placed  in  the  cylinder  edge  at  the 
previously  determined  crack  location  (see  figure  13).  Instituting  an  initial  crack  removes  the 
requirement  of  reaching  aps  to  fonn  a  crack.  The  energy  that  was  originally  applied  to  forming 
the  initial  crack  is  converted  to  crack  propagation  energy  instead.  Additionally,  to  exaggerate  the 
crack  growth  response,  the  Gc  values  are  reduced  to  1%  of  the  values  listed  in  table  1. 
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Initial  Crack 


With  these  modifications  included  in  the  model,  the  crack  propagates  farther  for  a  given 
deformation  and  achieves  complete  penetration.  Figure  14  shows  the  crack  length  after  a  plate 
displacement  of  0.125  inch  (3.18  mm),  matching  the  previous  analysis  and  experimental  test.  As 
the  crack  approaches  the  reinforcing  hoop  layer,  it  begins  to  turn  before  reverting  to  its  radial 
path.  This  crack  path  is  somewhat  similar  to  that  in  the  experimental  results,  where  the 
intralaminar  matrix  crack  turns  to  a  delamination  between  the  hoop  and  longitudinal  layer.  The 
delamination  fracture  toughness  is  lower  than  the  intralaminar  fracture  toughness  of  the  hoop 
layers.  In  this  model,  however,  the  difference  in  fracture  toughness  may  not  be  sufficient  to 
substantially  change  the  maximum  tangential  stress  direction  in  which  the  crack  propagates.  The 
crack  then  reverts  back  to  a  radial  crack  and  propagates  through  the  hoop  layers.  One  possible 
solution  to  this  modeling  delamination  as  a  separate  phenomenon  is  to  incorporate  a  cohesive 
strength  between  each  constitutive  layer,  which,  in  conjunction  with  the  strain  energy  release 
rates  available  in  this  analysis,  would  provide  an  additional  method  of  failure. 


Figure  14.  Failure  of  Cylinder  with  Strain  Energy  Release  Rates 
at  1%  of  Values  Listed  in  Table  1 
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CONCLUSION 


In  this  report,  the  extended  finite-element  method  (XFEM)  was  incorporated  into  a  3-D 
finite-element  model  to  simulate  lateral  compression  of  0°/90°  unidirectional  carbon  fiber/epoxy 
cylinders.  The  results  were  compared  with  failure  modes  observed  in  representative 
experimental  tests.  The  finite-element  model,  utilizing  XFEM,  replicated  the  test  parameters  and 
loading  procedure  by  instituting  a  prescribed  compression  between  two  flat  plates. 
Orientation-dependent  material  properties  were  used  for  each  of  the  10  distinct  composite  layers. 
Different  fracture  criteria  were  applied  to  the  longitudinal  and  hoop-oriented  layers  based  on  the 
respective  failure  mechanisms:  fiber  fracture  for  the  90°  hoop-oriented  layers  and  intralaminar 
matrix  fracture  for  the  0°  longitudinal-oriented  layers. 

Experimental  failure  initially  occurred  in  the  outermost  longitudinal  layer  of  the  cylinder 
and  propagated  through  the  matrix  radially  inward.  Upon  reaching  a  hoop-oriented  layer,  the 
crack  path  turned  to  a  delamination  between  the  hoop  and  longitudinal  layers.  In  the 
computational  model,  the  crack  fonned  in  the  outermost  layer  as  well  and  near  the  same  location 
as  that  of  the  experimental  test.  The  crack,  however,  did  not  fully  propagate  through  the 
composite  tube  as  was  seen  in  the  experimental  specimens.  An  initial  crack  was  instituted  into 
the  model  to  represent  surface  defects  in  the  composite  specimens.  Furthermore,  the  fracture 
parameters  were  decreased  to  facilitate  crack  growth.  When  these  changes  were  instituted,  the 
crack  path  more  closely  reflected  what  was  observed  in  the  experimental  specimens. 

This  research  demonstrated  the  ability  to  compare  failure  in  composite  materials  using 
XFEM  with  experimental  results — thus  improving  modeling  efficacy,  diminishing  the  need  for 
build-test-build  methodology,  and  potentially  reducing  cost  and  time  requirements.  Further  work 
must  be  done  to  accurately  account  for  material  defects  and  the  effect  of  localized  discontinuities 
on  the  fracture  toughness  of  unidirectional  composite  laminates.  The  goal  is  that,  with  more 
accurate  modeling  practices,  computational  analysis  of  composite  fracturing  can  reduce  the 
current  iterative  build-test-build  process. 
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